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Abstract 



We analyze an anisotropic spin 1/2 two legs ladder in the presence of 
various type of random perturbations. The generic phase diagram for the 
pure system, in a way similar to spin one chains, consists of four phases: an 
Antiferromagnet, a Haldane gap or a Singlet phase (depending on the sign 
of interchain coupling) and two XY phases designated by XYl and XY2. 
The effects of disorder depend crucially on whether it preserves XY symme- 
try (random field along z and random exchange) or not (random anisotropy 
or random XY fields). In all cases we computed the new phase diagram 
and the correlation length for the disordered system. The ladder exhibits 
a remarkable stability to disorder with XY symmetry. Not only the singlet 
phases but also the massless XYl phase are totally unaffected by disorder in 
stark contrast with the single chain case. Even in the presence of an external 
magnetic field breaking the spin gap most of the phase diagram (including 
the XY point) remains unaffected by the disorder, again in opposition with 
the single chain case. This stabilization towards disorder is similar, albeit 
stronger, to the one occuring for fermionic ladders. On the other hand the 
XY2 phase is very strongly suppressed. Disorder breaking XY symmetry has 
much stronger effects and most of the phases at the exception of the singlet 
one are now destroyed. Interestingly, the disordered XYl phase has a much 
shorter correlation length than the disordered XY2 phase contrarily with the 
case of perturbations preserving the rotation symmetry around the z axis. 
The ladder system is thus in fact much more anisotropic than its single chain 
counterpart with identical exchanges. Finally we examine the case of strong 
disorder or weak interchain exchange for XY symmetric disorder. Close to 
the isotropic point, when the interchain exchange is increased the system un- 
dergoes a transition between two decoupled disordered spin 1/2 chains to a 
singlet ladder unaffected by disorder. For more XY like systems, the transition 
occurs in two steps with a non disordered XYl phase occuring between the 
decoupled chains and the singlet phase. Comparison with other disordered 
gapped systems such as spin one chains is discussed. 
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I. INTRODUCTION 



One dimensional antiferromagnetic spin 1/2 systems are fascinating systems. They can 
present quasi-long range order, without symmetry breaking yet with critical spin spin cor- 
relation functions and diverging magnetic susceptibilitieJl. Not surprisingly, owing to the 
Jordan- Wigner transformation, such properties are remarkably similar to the ones of inter- 
acting one dimensional fermions and are known under the generic name of Luttinger liquid 
physics™. Experimental realization of spin 1/2 chain provide a good confirmation of these 
remarkable properties!'!. In addition, in a way reminiscent of one dimensional fermions , spin 
chains are extremely sensitive to disorder. In fact using the Jordan- Wigner transformation!, 
spin chains can be mapped onto one dimensional spinless fermions, where the coupling 
gives rise to an interaction term. Spin 1/2 systems are thus just another, and maybe more 
convenient experimental realization of disordered interacting fermions (or bosons)Bi. The 
transition to the disordered phase has been analyzed numerically in Refs. [T^-|T^ in good 
agreement with renormalization group treatments!'!. 

Recently it was realized that these remarkable properties get drastically modified when 
spin 1/2 chains are coupled together leading to the so called ladder systemEl. In this case, 
in a way very similar to the Haldane spin-S problem0'El, a gap is found to open for an even 
number of chains while the system remains massless for an odd number of chains. This 
phenomenon has been thoroughly investigated both analyticallyEUli and numericallyll3c3. 
In the recent years, progress in solid-state chemistry resulted in the apparition of com- 
pounds with a 2-legs ladder structure such as SrCu203, (La,Sr,Ca)i4Cu2404i and Cu2(l,4- 
Diazacycloheptane)2Cl4, in which the theoretical predictions about the gap have received 
an eclatant confirmationlih^l. 

Given the new physics present in these systems it is of paramount importance both from 
a theoretical point of view and for practical experimental systems, to investigate the effects 
of disorder in these systems. This is especially true for the two leg ladder system for which 
zinc doping experiments are already performed^. Due to the formation of a gapped singlet 
phase for the isotropic system it is somehow natural, although not always trivial, to expect 
that weak disorder will have no effect. Most of the studies about the influence of disorder 
have been thus devoted so far to the effects of strong disorder on such a gapped phase either 
in ladders or in other gapped systemsS^ However, in the spin 1 chain system, the gap is 
proportional to the intra chain exchange, whereas in the ladder system, the gap is directly 
controlled by the interchain hopping J^, and thus can be made small at will compared 
to the intra chain exchange J permitting the use of perturbative methods to analyze the 
competition of randomness and singlet formation. Moreover, very little is known on the 
effects of disorder on the anisotropic system, and inparticular in its gapless phases. A 
surprising effect found for fermionic two lee laddera3E£l in which a gap opens in all the 
modes except the symmetric charge modeE3E3 is a remarkable stability of the ladder system 
to disorder. Whether a similar stabilization by gap opening exists for two coupled spin 
chains needs to be investigated. We undertake such a study in the present paper. Quite 
remarkably we find that an anisotropic spin ladder is completely stable even in the gapless 
phase any weak disorder respecting the XY symmetry. This is in stark contrast with a single 
chain or even with ladders of spinless fermions. 

The plan of the paper is the following: In section |I|, we recall the basics of the bosoniza- 
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tion treatment of spin chain systems. In order to study the effects of disorder on the 
anisotropic spin ladder we first investigate in details in Section |T| the phase diagram of the 
pure system. In addition to the three phases discussed in Ref. |16|, namely a Singlet phase, 
an XYl phase and the antiferromagnetic phase, we find another XY phase, the XY2 phase, 
in analogy with the case of spin one chains^. We discuss the possibility of observing this 
XY2 phase at intermediate coupling. In section |I^ we investigate the response of the two 
chain system with respect to weak disorder preserving the rotational symmetry around the 
z axis, i.e. random fields parallel to the z axis and random exchange. The spin ladder shows 
a remarkable stability to these perturbations even in the gapless phase. We also consider 
the disordered ladder under a magnetic field, able to suppress the gapped phase. Finally we 
study the crossover between a single chain and the ladder behavior obtained by increasing 
the interchain coupling. Given the remarkable stability of the ladder to random perturba- 
tions respecting the XY symmetry it is interesting to also investigate that breaks rotational 
symmetry around the z axis. We thus study in section |V| random anisotropy and random 
fields in the plane perpendicular to the z axis. Quite surprisingly the effect of these pertur- 
bations is opposite to the one respecting the XY symmetry: the XYl gets very unstable, 
whereas the XY2 phase is much less affected. Finally conclusions and open questions can 
be found in Section |V|. 



II. BOSONIZATION OF ONE DIMENSIONAL SPIN 1/2 CHAINS 

In this section, the bosonization of spin 1/2 chains is briefly recalled. For the sake 
of definiteness, an XXZ spin chain will be considered. It is described by the following 
Hamiltonian: 

^ i i 

Spin chains Hamiltonians can be transformed into interacting 1 dimensional fermion systems 
by expressing the spin operators S'+,S'~,S'^ in terms of fermion operators a\a using the 
Jordan Wiener transformationHi'il 



Sf = (— )*aj cos n ^ a^aj 



j=0 

Sj^ = (— )* cos TT ^ Ojaj ai (2) 
V j=o J 

Si = alai - 2 

This changes the XXZ model into a model of spinless fermions with nearest neighbor inter- 
action described by the Hamiltonian 

H = -Jx/2^(4a„+i + ai^^ttn) + Jz^{4(^t - - ^) (3) 
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Jz = corresponds to free fermions and the zero magnetization sector corresponds to an 
half-filled band with Fermi- wavevector kp = a being the lattice spacing. 

To treat the problem with a finite J^, it is convenient to use the boson representation of 
one-dimensional fermion operators. Only the final results are eiven here since details about 
the method can be found in various places in the literaturei^@. A detailed derivation with 



the same notations can be found in the appendix A of Ref. |36|. The spin operators are given 
by 



Sz{x) 



-ie{x) 



V2 



na 



e + cos20(a;) 
cos2(j){x) 



TT na 



(4) 



Where S'^(x) = S^{x) = ^ for x = na, a being the distance between two nearest 
neighbors sites along the chain. The phase is related to the average density of fermions 
(or equivalently to the uniform spin density along z) by Sz{x) = —^dx(p, whereas 9 is 
connected to the conjugate momentum 11 of (such that [0(a;), n(a;')] = i6{x — x')) by 
9{x) = /j^oo dyll{y)- In a very crude sense (f), 9 can be viewed as the polar angles of a spin. 
The Hamiltonian (|l]) gives 



xxz 



dx 
2^ 



+ 



2A 



(27ra) 



dx cos(4(/)) 



(5) 



This is a sine-Gordon Hamiltonian with /? = 4. The cosine term comes from the so 
called Umklapp process. Such processes are possible as ikp is equal to a reciprocal lattice 
vectoiiM. The coupling constant can be determined in perturbation in the interaction 
Jz and gives for Jz -C Jxy 



K 



1 + 



4X 



-1/2 



TT J, 



u = aJx 1 + 



A = —Jztt 



4J. ^ 

7^ Jxy, 



1/2 



(6) 



The range of validity of (^ is much larger than the simple perturbative regime provided the 
correct constant K and u are use dW. The isotropic point Jz = Jxy (Heisenberg model) 
corresponds in the bosonization description to K* = h,A* = as can be seen from the 
scaling dimensions of and 5*^ (see Eqs. ( |A^ ) and (|A3D of appendix 0). These correlation 
function become identical, as they should for the Heisenberg point, only for K = ^. At 
that point cos(4(/)) is marginally irrelevant. For larger values of Jz {Jz > Jxy) the cosine 
term becomes relevant. A gap opens in the excitation spectrum of the (p field, leading to 
an insulating charge density wave state for the fermions and equivalently to an Ising order 
along z for the spin chain (see e.g. Ref. p9|-|i2|). 

The 



Another instability of the XXZ model can be guessed from (|g). For Jz — 
K —>■ oo. This is the sign of a transition to a fully polarized ferromagnetic state. 
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exact solution shows that the ferromagnetic transition takes place for Jz = —Jx^^- In the 
following, ferromagnetic systems will not be considered as their behavior is very similar to 
the one of classical systems at T = QK. 



III. TWO COUPLED SPIN 1/2 CHAINS 

In this section, two XXZ chains coupled by an exchange term are considered using 
bosonization techniques. The exchange coupling is of the form : 



-n interchain — / , ^ 



(7) 



i,a=x,y,z 



where for simplicity, as for the intrachain Hamiltonian, planar isotropy J± = J± = J±^ is 
assumed. 

The total Hamiltonian is : impure = Hxxz,i + Hxxz,2 + -f^interchain- Each single chain 
Hamiltonian can be expressed in terms of fermions operators using the Jordan- Wigner trans- 
formation like in section |T|. One has to take care that the "fermions" operators should have 
anticommutation relations inside one chain but commutation relations between the chains 
to respect the spin commutation relations. This causes no modification of the bosonized 
expressions @ for the spin operators. Indeed it is easy to check that (^) preserves all the 
correct spin commutation relations. Hxxz has been obtained in section |l| and only the 
bosonized expression of -f/iaterchain is needed. Using (H) and keeping only the most relevant 
operators one obtain JlUllcJ: 



interchain 



2gi 



{2na) 



■ COS( t^i 



-Jfa / dx 



'292 



(27ra)' 



■ cos 2( 



+ 



2^73 



{2na) 



cos 2(01 + 



dx 



(8) 



Where : 



92 
93 



J la 
Jfa 



(9) 



The total Hamiltonian ifpure is rewritten in terms of the fields 



givmg 



^and 



h+4>2 
V2 



Hpuvc Hg -\- Ha 



7/ 



+ 



292 



(27ra)2 

293 
(27ra)2 



dx cos{V8(j)s) 
dx cos 



291 



{2na) 



(10) 

dx cos{V26a) 



Where : 
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/ KJla\ KJla 
Us = u{l^ — , Ks = K [1 — 

u^K being given bv 10) for a XXZ model and small J^, and more complicated expressions 
in the general casecU'ES. If the correct values of u and K are used (0) is valid irrespective 
of the values of Jz and J^y and only assumes J^a -C Because of the choice of an 
hermitian string operator, the cos(40i), cos(402) terms in iJxxz,i and -ffxxz,2 respectively 
can be dropped, since they always give terms that are less relevant than the terms that come 
from -^interchain- Notc that the interchain coupling in ( [I0| ) is different from the one existing 
in spinless fermions ladders0'il'il. This is due to the presence of the string operator 
in the fermion-spin equivalence. Therefore two coupled spin chains have a different physics 
that cannot be directly deduced from the fermionic resultSl. This distinction between 
spin chains and spinless fermions will be crucial in the presence of disorder. 

Eqs. (P) and (p!OD have been derived, for particular cases or for physically related sys- 
tems, in various places in theliterature. Such Hamiltonians were first introduced to describe 
a quantum spin-one chainli3S. For this problem, although the form of the Hamiltonian are 
identical g2i gs are related to parameters with a physical meaning different of the one 
of J±. For two coupled spin chains and (|TD|) have been obtained and invest igated0@ 
for the particular case of isotropic interchain coupling (J^'^ = J]_). This corresponds to 
taking g2 = g^ = gi/n in (p!oD . More recently, another derivation for isotropically coupled 
spin chains has been obtained by a different methodEl. Instead of using a Jordan- Wigner 



transformation, Ref. |18| uses a representation in term of localized fermion. Such a represen- 
tation is particularly well adapted to describe an isotropic coupling since the superexchange 
leading to the spin Hamiltonian is automatically generated. Interestingly one recovers the 
same bosonized Hamiltonian than (|I^) (with = 1/2 due to isotropy), but with a different 
relation between the coefficients g. In Ref. |18| the coefficient satisfy g2 = gs = gi J±. 
This difference comes from the fact that in the present derivation the calculation of the 
coefficients g is perturbative and the expression (|]) is accurate for <C J^y. In going to 
the isotropic point these coefficients get renormalized by irrelevant terms that do no affect 
the form of the low energy Hamiltonian (|10[), but can change the explicit value of the coef- 
ficients compared to @. Comparison of the two limits Jz -C Jxy and Jz ~ Jxy shows that 
this renormalization effects are rather weak, and that (|TO|-|ri|) provide an extremely accurate 
description of the physics of the ladder system for all values of Jz, provided one does not 
need the precise values of the coefficients g. 

Using ( p!OD the phase diagram of the pure system can be derived. Since the main goal of 
this paper is to analyze the additional effects of disorder, we restrict ourselves for the study 
of the pure system to isotropic interchain coupling but the intrachain coupling is arbitrary. 
Nevertheless, (0) can be used in more general cases. The derivation is technically similar 
to the one for the spin one chainEl, but gives of course different physical phases. Since Hg is 
a standard sine-Gordon Hamiltonian (see appendix 0) its spectrum has a gap if i^'^ < 1 and 
is gapless otherwise. The treatment of Ha is a bit more subtle. cos(A/80a) and cos{y/26a) 
have respective scaling dimension of 2Ka from ( [A3|) and {2Ka)^^ from (|A4|) so that both 
are relevant for 1/4 < Ka < 1. However, the coefficients of these operators are of the same 
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order of magnitude. Therefore, the most relevant operator is the first to attain the strong 
couphng regime under renormahzation. In the strong couphng regime this operator takes 
a mean value that minimizes the ground state energy. 9a and (j)a being conjugated, when 
one of them develops a mean value the exponentials of the other one have zero expectation 
values and exponentially decaying correlations (see appendix]^). Therefore, (pa acquires a 
mean value while 9a is disordered for Ka < 1/2 and the situation is reversed for Ka > 1/2. 
As a consequence four different phases exists, similarly to the spin one case@, regardless of 
the sign of the interchain coupling. Three of these four phases correspond to the one derived 
in Ref. but, as discussed in detail below an additional massless phase can exist. For 

isotropic intrachain coupling the results of Ref. |T8| are of course recovered. A summary of 
the phase diagram is shown on Figure [l|. Let us now discuss each phase in details. The 
physics of these phases depends whether interchain coupling is ferro or antiferromagnetic, 
since this changes the average value taken by the massive fields. 



A. Ferromagnetic interchain coupling 

A summary of the physics of each phase can be found in Table |. A qualitative under- 
standing of the results of table | can be obtained in the limit | |^ J^Jz- Then, the 
ground state energy is minimized by forming a spin one on each rung with the 2 spins 1/2. 
That case is thus physically identical to the spin one antiferromagnetEl. The same physics 
hold even if J±_ <^ J, Jz- As shown in Table |I| four sectors exist. 

In Sector I, the two spins 1/2 on every rung point in the same direction parallel to the 
z axis while each chain is antiferromagnetically ordered giving an effective spin one Ising 
Ant if er romagnet . 

In sector II a completely gapped phase exists where all spin-spin correlation functions 
decay exponentially. Using the analogy with the antiferromagnetic spin 1 chain, that phase 
is identified with a Haldane gap phase0. Singlets are formed here along the legs of the 
ladder^. The identification can be made more complete by proving that the VBS order 
parameter, characteristic of a Haldane phase, is effectively non-zero (see appendix 0). 

Sector IV is a phase where only XY correlations remains. This XYl phase has order 
parameter Si + S2 = in agreement with the physical picture of a total spin equal to one. 
Semiclassically, the XYl phase tends to have the two spins in the XY plane. This phase 
has only a gap in the antisymmetric spin excitations and has gapless symmetric excitations. 
(See table |, sector IV for the bosonized expression of the order parameter). 

Finally, in sector III, S1S2 which can be viewed as an effective raising operator for a 
spin 1/2 is an order parameter. Its bosonized expression can be found in table |. In the case 
of a true spin 1 chainEl, such a phase can be obtained by adding an anisotropy of the form 
—D{Sf + 5*1)^ which favors the states = +1, —1 over = 0. Thus, the XY2 phase is a 
phase in which the two spins on each rung can only be found with either {Si, S2) = (1/2, 1/2) 
or [Si, S2) = (—1/2, —1/2). They can make transitions from one of these state to the other 
and therefore form an effective spin 1/2. This accounts for the critical fluctuations in that 
phase. The XY2 phase also shows subdominant Ising Antiferromagnet fluctuations due to 
the existence of the intrachain antiferromagnetic coupling. All other single spin fluctuations 
decay exponentially fast. 
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At small J_L, (0) gives Ks ~ Ka- Thus, only the Haldane gap, XYl and Ising An- 
tiferromagnet phase can be reached. These three phases were studied in Ref. 16,44. For 
isotropic coupling refermionization indeed gives an Haldane phaselll. When is very large 
the effective Hamiltonian for the resulting spin one chain contains only J^y, Jz terms, and 
thus only the same phases can be reached. However, since when increases, (^Tj) indicate 
that Ka decreases, whereas Ks increases it is interesting to investigate whether one could 
reach the XY2 phase an intermediate J±. Of course answering this question probably needs 
a numerical study of the model, since all the couplings would be of the same order. Even if 
this interesting possibility does not occur, the XY2 phase can always be obtained by using 
more complicated interchain couplings. 

Finally, it is noteworthy that the only phase with SU(2) symmetry is the Haldane gap 
phase. Two ferromagnetically coupled isotropic spin 1/2 chains, thus always form a Haldane 
gap phase, irrespective of the strength of the coupling. The other phases are characteristic 
of anisotropic systems. 



B. Antiferromagnetic interchain coupling 

Results for the antiferromagnetic coupling are summarized in table |T|. 

Sector I is again an Ising antiferromagnet phase. The spins on the same rung point in 
opposite directions so that in lattice representation, the order parameter is {—Y{S^ i~Sn2)- 
The boson representation of this order parameter can be found in table 

Sector IV corresponds to the XYl phase obtained in Ref. |l^. Its order parameter is 
Si — S2 ■ See table |1| for the bosonized expression. In that phase, the spins of each rung 
stay in the XY plane and are antiparallel. They are free to rotate around the z-axis provided 
they remain antiparallel and are antiferromagnetically coupled to their neighbors giving rise 
to an XY phase. 

Sector HI is the other XY phase, the XY2 phase, for which the order parameter is 
SiS2- The hidden spin 1/2 degree of freedom corresponds now to the two states of zero z 
component of the spin formed by the addition of the spins 1/2 on the same rung. Tothe 
best of our our knowledge, such a phase was not investigated in the previous studiesei'lii. 
The XY2 phase presents subdominant Ising Antiferromagnet fluctuations. Thus, in the XY2 
plane, there is a tendency to have the 2 spins parallel to the z axis and antiparallel among 
themselves. The intrachain coupling allows transitions between these two states leading to 
critical fluctuations analogous to the one of a spin 1/2 system. 

Finally, there is a Singlet phase in which there is a gap to all excitations and all correlation 
functions of spin operators decay exponentially fast. A simple picture of that phase is 
obtained in the limit J±/Jz — ^ 00. In that limit, the ground state can be seen as made 
of singlets along the rungs. In such a phase, there is an effective zero spin on each site 
and accordingly all the spin-spin correlation function are zero. For finite J±/Jz, massive 
triplet excitations can propagate along the chain leading to exponentially decaying spin- 
spin correlation functions. 

As in the case of a ferromagnetic interchain coupling, at small J±/Jz the only possible 
phases are the XYl phase, the singlet phase and the Ising Antiferromagnet phase analyzed 
in Ref. |16|. The XY2 phase appears at Kg > 1 and Ka < 1/2 (see table H) and thus. 
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as for the ferromagnetic case, cannot be obtained for purely local intrachain coupling and 
small Jj^. Moreover, for antiferromagnetic interchain coupling, increasing does not favor 
the formation of the XY2 phase, so it is likely that one needs here longer range interchain 
coupling. Finally, it is important to remark that as in the case of the ferromagnetic interchain 
coupling, the only phase that has spin rotational symmetry is the singlet phase, and that 
all other phases can only exist in systems with anisotropic interactions. 

IV. RANDOM FIELDS AND RANDOM PLANAR EXCHANGE 

Since we now have a complete description of the phases in the pure system, the effects of 
various type of disorder on these phases can be considered . Very different physical effects 
occur depending on whether the disorder respects or breaks the U{1) symmetry (i. e. the 
rotational symmetry around the 2;-axis) of the system). For the case of symmetric couplings 
in the XY plane for the pure system the first type of disorder is the one most likely to occur 
physically. We consider its effects in the present section. As we discuss in more detail later 
perturbations that break the U{1) symmetry are nevertheless also interesting to study and 
are considered in section 0. 

A. Coupling to disorder 

Two type of U{1) conserving disorder can be investigated. The first one is the standard 
random field along the z direction, while the other is the random exchange. If a representa- 
tion of the spins in terms of hard core bosons was used such perturbations would correspond 
to random potentials or random hopping along the chains for the associated disordered boson 
problem^. 

1. Random field along z 

Let us consider first the coupling to a random magnetic field along the z direction. 
This field is assumed weak enough not to destroy the gaps of the pure system. The lattice 
Hamiltonian for that problem is if = impure + Hzf with 

HzF = [^A^li + K^^l^ 

i 

For simplicity the disorder is taken Gaussian and uncorrelated from chain to chain and 
from site to site: 

h!Jhl2 = (13) 
hfj^ = Dkj (14) 

First, a bosonized representation of (|12|) has to be obtained . To do so, one goes to the 
continuum limit, decompose the fields in a slowly varying {q ~ 0) and an oscillating part 
(g ~ 2kp = ^), introduce the bosonized expression for spin operators (§), and keep only in 



(12) 
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the coupling to disorder the slowly varying terms, in a way similar to Ref. |],^. The final 
result in terms of the 0i and 02 fields is : 



H, 



ZF 



dx [/if ^(x) cos20i(x) + /if ^(x) cos202(x) 



(15) 



with 



/ij ''{x)hj -^(x') = D2kpa'S{x — x')6ij 
hr\x)h]^\x') = Dg^oaS{x - x')6ij 



(16) 



All other correlators are zero. Introducing the symmetric and antisymmetric combinations 
/i^~'^(x) = hf"^{x) + /i2~°(x), and h'^^^{x) = hl""^{x) — h1"'^{x), ([T5|) is rewritten in terms of 



the fields 0^, 0^ as : 



h: 



ZF 



ZF 



-"ZF "T -"ZF 



dx 



TT 



/ir°(a;)40s(x) + /ir°(a;)9,0„(x) 



-"ZF 



/if ^ (x) cos v^(0a + J (x) + +/if ^ (x) COS v^(0a " 0s) (a^) 



(17) 
(18) 

(19) 



(|T^) is the starting Hamiltonian for the discussion of the effect of a random z field in the 
framework of bosonization. ( [I8| ) represents the coupling of the system with the g ~ 
components of the random field, and (|19]) the coupling with the 2kp components. 

The dxCps term in ( [T^ ) has to be discussed first. Two cases have to be distinguished : 
when there is no gap in 0^, a transformationll 0s — > 0s + dx'hf"^{x'), results in a new 
Hamiltonian of the form ( [T7| ) with /i^~° = 0, and identical correlators for the disorder. 
Therefore, the hf"^ terms cannot affect the low energy spectrum of the system and only 
affect its correlation functions in a trivial way. On the other hand, if there is a gap in the 
0s degrees of freedom, through refermionization the problem in the 0s sector is equivalent 
to a filled band of fermions with a gap to excitations in the upper band and a small slowly 
varying random chemical potential. The random potential being smaller than the gap, it 
cannot induce transitions in the upper band (i.e. the system is incompressible) and thus 
cannot disorder the system. Therefore, in both cases, at weak disorder strength the /i^~° is 



always irrelevant. A similar situation occurs for the dx4>a term in ([T8|). When 0^ has a gap 
the problem is obviously identical to the problem of the gap in 0s. When 9a is massive, the 
refermionization procedure shows that computing the response to the dx4>a term amount to 
computing the superconducting response function of a band insulator. Clearly, this has an 
exponential decay so that the random field coupled to dx4>a term cannot disorder the system. 
Thus all the g ~ terms (^8[) have no important effect on the low energy response of the 2 
chain system and can be discarded. If the disorder is stronger however, the g ~ part of 
the disorder will however be able to destroy the gaps of the pure system. This case will be 
discussed later. 

Having shown that the g ~ do not affect the phase diagram at weak disorder, the 
2kF terms (ITPI) remain to be treated . To simplify the notations the shorthand hi is used 
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for hi'''' . The effect of the 2k p disorder depends on the phase that would exist in the pure 
system in the absence of disorder. A detailed phase by phase discussion of the effects of 2k p 
randomness is given in section |IVB| . The phase diagram in presence of disorder is shown in 
figure 0. 



2. Random exchange 



The other type of disorder that preserves the rotational invariance around the z axis is 
random planar exchange and random z exchange. For convenience the random z exchange 
and random planar exchange are treated separately. 

Restricting to the small disorder case, the Hamiltonian in the planar exchange case is 

H = impure + -f^PE with 



TT _ \ qx j_ qy qy \ j_ r'^f q^ qx , qy qy 

-npE — 2^ [-'n Wn,l^n+l,l ^n.l'^n+l.lj T "^n Wn,2^n+1,2 ^n,2'-'n+l,: 



(20) 



with JnJn' = D6n,n'Sp^p^ p = 1,2 being the chain index and <^ J. The bosonized 
Hamiltonian is given bj 



PE 



dx 



Ji{xf^^ sin 201 (x) + J2{xf^^ sin 202 (a;) 



+ / dx^^Uxf-' 



(21) 



where Jp{x) is the continuum limit of J^. 

Let us now consider the random z exchange case. In the lattice spin 1/2 representation, 
the coupling to disorder is represented by 



random Jz 



El q^ q^ _L 7 q^ q^ 



(22) 



The bosonized form for the g ~ term is easily obtained in the form 



(g~0) 
random 



dx 



J': 



-0,1, 



x){dAinx) + jr'^\x){dM\^ 



(23) 



Obtaining the 2kF is not so straightforward. The derivation is given in Appendix 0. The 
final result is 



random 



dx 



j2fc^,i(-^^) Si^20i(x) ^ j2kp,2, ^sin202(x) 



27ra 



2'Ka 



(24) 



with Jf^p{x)J'^^p{x') = D5{x — x'). The breakdown of bosonization for too large random 
planar exchange can be read off in the g ~ terms: if g ~ disorder becomes too large 
the quadratic part of the Hamiltonian is no more positive definite leading to a breakdown 
of the bosonization description. However, if the weak disorder condition is met the g ~ 
terms couple to H^ and (9x0)^ terms. Power counting then implies that the g ~ terms are 
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irrelevant. Thus, the stabihty against weak random exchange is again determined only by 
the 2k F terms. 

Bosonization leads to the same Hamiltonians ( pT]) and for the coupling with the 2k p 



part of random planar exchange and random z exchange. These two perturbations can thus 
be treated the same way. In addition, in the phases in which 0^ is gapless, one can make 
the transformation: 



^2 



02 + vr/4 (25) 



changing the random exchange term into a random z field as can be seen by comparing 



and (pTD for the random exchanges and (19) for the random field. Thus in the XYl and 



XY2 phases both perturbations lead to the same stability regions and the same correlation 
lengths as is discussed in Section |1V B| . Of course the physical properties of the disordered 



phase are different depending whether the disorder is random exchange or random field. The 
phase diagram in the presence of random exchange is given in Figure H. 



B. Effect of disorder 

In this section the effects of both the random 2;-field and the random exchange are 
considered in detail . The simpler case of the XY phases is discussed first, then the case of 
the Haldane gap phase and finally the case of the Antiferromagnet phase. 



1. XY2 phase 

The XY2 phase i. e. The X-Y phase with non-zero mean value of (pa (sector III of 
tables I and H) is considered. For the sake of definiteness, only the random z field case 
is analyzed , but as discussed before in the XY phases the same results also apply in the 
random exchange case. Using (|1^) the effective Hamiltonian describing the coupling with 
the impurity potential is 

HzF = C J dx{hi{x) + h2ix)) cos(y20,)(x) for < (26) 
= C' J dx{hi{x) + h2{x)) sm{V2(f),){x) for > (27) 

where C = (cos 0a) and C = {sin (pa) (see table | and table The relevance of disorder 
can be determined by looking at the renormalization of the disorder termi. 

The RG equation for disorder is obtained from the scaling dimensions of operators en- 
tering the coupling with disorder as it is the same for both ( p6D and (|27|) . 

^^ = {3-K,)D'''-{l) (28) 

Note that the RG equation for the disorder does not depend on the nonuniversal constants 
C, C. The RG equation for the Luttinger liquid parameter Ks is : 



12 



(29) 



where A = C, C" depending on (0a)- This equation is non-universaL However, for very 
weak disorder the region of stabihty is given by (^Sj) , and one can discard (|29|) except to 



compute the critical properties very close to the transition^. Even for finite disorder 
gives correctly the renormalized value of the parameter K, and thus the exponents of the 
correlation functions at the transition. 



Equation (pHj) immediately shows that the XY2 phase is unstable unless Kg > ?>. It 
also gives the correlation length in the disordered phase, for weak disorder and if one is not 
too close to the transition point!. The renormalization of Kg at weak disorder, (pQ]) being 
negligible, the RG equation for D can be integrated into D{1) = D(0)e^^~^*-". This form of 
D{1) is valid as long as D <^ —, since this condition ensures that the energy scale induced 

by disorder is much smaller than the energy cutoff. For D{1) ~ —(0) disorder cannot be 
treated as a perturbation and the RG has to be stopped . At that lengthscale, the system 
appears to be strongly disordered. Under the RG flow, the cut-off length has increased from 
a(0) to a{l) = e'a(O) at which point the system appears to be strongly disordered . Below 
it the system is described by bosonization and disorder can be treated as a perturbation. 
Thus the length a{l) is the correlation length ^, leading to 
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Using the transformation (pSD , the same region of instability Kg < 3 and the same correlation 
length ^ = (l/D)3-ifs are obtained for the random exchange case. 



2. XYl phase 

The case of the XYl phase in the presence of a random z field is now considered . This 
case is more involved since to first order the effective coupling with disorder is zero (see 
tables H, |, sector IV). However, the coupling to disorder although irrelevant by itself can 
generate relevant terms in higher order in perturbation. Such terms can be computed, as was 
done in the case of the two chain of fermions problemUS, by integrating over the massive 
modes. This leads to the effective Hamiltonian for the coupling of the massless modes 0s to 
disorder 

-f^zF.cff. = j dxrjcsix) cos{\/8(j)s{x)) (31) 

with ries{x)ries{x') oc D'^6{x — x'). By using the same method than above, the scaling 
dimension of the disorder is now 3 — 47^^, implying that the disorder is irrelevant unless 
Kg < ^. But a.t Kg = 1 the g2 term becomes relevant and drives the system towards the 
singlet (if g2 > 0) or the Haldane gap (if g2 < 0) phase (see tables |T|J3). Thus, the XYl 
phase is unaffected by weak random 2;-fields (or using the transformation (^) by random 
exchange) since these perturbations are irrelevant in its whole domain of existence. The 
XYl phase is thus unaffected by all the perturbations respecting the rotation symmetry 
around the z axis. 
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3. The gapped phases 



From section |1V B 2| one sees easily that a two leg ladder system in a singlet or a Haldane 
gap phase coupled to random z fields is described by the following effective Hamiltonian 
coming from integration of the massive antisymmetric modes 



H= ^ 



dx 
2^ 



+ dx 



2g2 + VcsXx) 
(27ra)2 



cos V 8 



(32) 



In the above discussion of the XYl phase, it was shown that the non-random g2 term is more 
relevant than the disorder term. Thus if the disorder is weak enough the deterministic part 
of the cos a/80s term dominates over the random one and the gapped phases are stable with 
respect to small random z-fields or small random exchange. Such stability is reasonable on 
physical grounds in a singlet phase. The stability to stronger disorder will be discussed in 



Section IV C 3 



4- The Ising Antiferromagnet 

The Ising Antiferromagnet phase (sector I in tables |,||) is the simplest to discuss due 
to its classical character. However, a marked difference appears depending on whether one 
considers a random z field or a random exchange. For J_|_ < one has {{SI + S'O) ^ and 
for J_|_ > 0, {{SI — Si)) 7^ 0. For the random z-field, a simple Imry-Ma type argument^ 
shows that the long range order is lost and that the ground state of the disordered 2 chain 
system is made of domains of characteristic size ~ This result is identical to the 

case of a classical antiferromagnet. It can also be obtained directly on the spin Hamiltonian 
without using bosonization. The Ising Antiferromagnet phase of the two leg ladder system 
is thus unstable in the presence of an arbitrarily weak random z field, as was the case in the 
one chain problemiS. On the other hand, in the case of a weak random exchange, there is 
no coupling at all to the disorder and the Ising Antiferromagnet is stable in the presence of 
disorder. 

A summary of the phase diagrams for the random z field and random exchange are 
respectively shown in Figure ^ and Figure ^. 



C. Discussion and Physical properties 

1. Comparison with a single chain 

It is interesting to compare the above results with a single disordered spin chain. For a 
single chain only two phases exist in the pure system: the Ising antiferromagnet (for K < 
1/2) and an XY phase (for K > 1/2) that is the analogous of the XYl phase for the ladder 
system. The XY phase is destroyed by a random field along z or random exchange for K < 
3/2ii'0. The antiferromagnetic phase is unstable in the presence of a random magnetic field 
but not random exchange as shown by an Imry-Ma argumentiS. The isotropic Heisenberg 
point is thus unstable to infinitesimal disorder. 
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For the ladder case the isotropic point corresponds now to a singlet gapped phase. Quite 
naturally this phase is insensitive to small disorder. The antiferromagnetic phase gives 
identical results than for a single chain. Quite remarkable however the XYl massless 
phase is now resistant to all perturbations respecting around the z axis, at the opposite of 
the corresponding phase for the single chain. This surprising result can be understood by 
noticing that the XYl phase is much more anisotropic than its one chain counterpart (in 
particular although it is massless it still has exponentially decaying 2kp correlations of the 
Sz component). The disorder can only couple to higher operators that are less relevant, in 
a way reminiscent of the situation in the fermion ladder problemilS. 

The opposite situation occurs for the XY2 phase. This phase phase reveals itself much 
more unstable than the XY phase in the presence of random z fields or random exchange. Its 

domain of existence is reduced to Kg > 3. As a result, it correlation length ~ [j^ ^'^^ 

is also much shorter than its single chain counterpart l\ ch. ~ [j^ ^"^^ • This is due to the 
existence of strong antiferromagnetic fluctuations in the XY2 phase that easily couple to 
disorder. It is again reminiscent of the easy coupling of the charge density wave phase for the 
fermionic ladderBH. The two most remarkable effects occurring in the ladder system are 
thus the stabilization of the isotropic point due to the singlet phase and quite unexpectedly 
the stability of the XYl massless phase as well. This remarkable stability against disorder 
prompt for several questions. In particular it is interesting to understand how the phase 
diagram evolves if the gaps allowing for this stability against disorder are destroyed, either 
by a magnetic field or by increasingly strong disorder. These two cases are examined in the 
next two sections. 



2. Effect of a magnetic field 
The coupling to a uniform magnetic field has the form: 



-f^uniform = / dxd^(f)s (33) 



TT 



Therefore, a strong enough magnetic field, suppresses the gap formation in (ps but does not 
affect the antisymmetric sector. Thus the antiferromagnetic phase and the Haldane gap 
phase disappear upon application of a strong enough magnetic field and are replaced by 
incommensurate phasesea. The stability of the remaining two incommensurate gapless XY 
phases in the presence of the random magnetic field or the random exchange can be analyzed 
by similar techniques than the ones of section [IV B. 



Using the analysis of Sec. [IV B 2| , one finds that the random z field results in the sup- 



pression of the XYl phase for Kg < 3/4 and formation of a random antiferromagnet with a 

2 

correlation length lAF^ikp ~ {j^ 'i-iKs ^ rpj^^ modulation of the spin density is at a wavevector 
Akp = 27r(l + 2m) , where m is the magnetization and is incommensurate^ with the the 
lattice spacing. Contrarily to the zero magnetic field case, the disordered Akp antiferromag- 
net is no more wiped out by Haldane gap formation. According to Sec. [IV B 1[, a random 



z field also results in suppression of the XY2 phase for Kg < 3 and formation of a random 
2kF antiferromagnet . The situation in this case is however identical to the zero magnetic 
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field case. The results are summarized on figure |^. In particular, for 1/2 < Kg < 3/4 the 
system is again sensitive to disorder. The case of a random exchange is similar except that 
the Akp antiferromagnet is replaced by a Akp random singlets phase. 

Using the standard mapping of spins on hard core bosons, the problem of the coupled 
chains under field can be related to the more general problem of a bosonic ladder. This 
problem will be analyzed in details elsewherecj. In the boson language, the XY phases 
correspond to superfluid phases, and antiferromagnetic phases to charge density waves. 
Random z fields correspond to random potentials and random planar exchange to random 
hopping. 



3. Stability to strong disorder 

It is also interesting to look on how the gaps in the system can be destroyed by the 
disorder itself, when it becomes strong enough. Strictly speaking such a study is beyond 
the reach of a renormalization group treatment since such a transition would correspond 
to a transition between two strong coupling fixed points , in a regime where the RG is not 
rigorously applicable. Fortunately one can still make some physical arguments. It would be 
thus of great practical interest to check whether the simple analysis performed here can be 
confirmed in more sophisticated treatments such as simulations or by using non-perturbative 
method&i'ffl. 

For simplicity we concentrate on the Haldane (or singlet) phases and on the XYl phase 
in the presence of a random magnetic field. Both the symmetric and antisymmetric gap, 
when they exist can be destroyed by either the g ~ or g ~ 2kp component of disorder. 
A similar effect has been analyzed in the case of non-magnetic impurities in Ref. |5^ using 
bosonization and RG techniques. The g ~ component of the magnetic field can prevent 
(pa from developing a gap, however if 6a develops a gap, (i. e. Ka > 1/2) this component of 



the random magnetic field has simply no effect for the reasons exposed in Sec. |1V A . 

Three regimes of stability can thus a priori be defined. At very small one has disor- 
dered decoupled spin 1/2 chains where no gap exists both in the symmetric and antisym- 
metric sector. The correlation length in this regime is ^decoupled ~ (l/-D)^/^^~^'^\ using the 
approximation ~ Ka ~ K for small J±. This phase is unstable when the correlation 
length due to the opening of one of the gaps by J±_ becomes comparable to ^decoupled- It is 
easy to check that the shortest correlation length induced by is the antisymmetric one. 
The transition occurs therefore when 

J* ^ ^(2-l/2i^)/(3-2X) ^3^^ 

thus giving J* ^1/2 for K = 1/2 (the isotropic point) and Jl ~ D^'"^ for i^' = 1 (the limit 
to the XYl phase). For larger values of J±_ a gap exists in the antisymmetric 6a mode. 

Since the gap in the symmetric mode behaves as jY'"'^ the properties of this mode 
crucially depends on the value of K. The critical value of J± needed to resist the 2kp and 
the g ~ of the random z field of the disorder are respectively 

J±,ci ~ Dlzr^ (35) 

i-iKs 

J±,c2 ~ D'^ZZ (36) 
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For the isotropic point K = 1/2 both these values are smaller than J]_ . By increasing one 
should thus go directly from two disordered decoupled spin 1/2 chains to a stable Haldane 
(or singlet) ladder system. Although one can only give physical arguments for the nature of 
the transition, the crossing of two correlation length suggest that there is a sudden drop of 
the gap to zero. It is clear that the decoupled chains do not have any form of topological 
hidden order and that their response function is simply the one of isolated disordered chains. 
The phase diagram is shown on Figure ^a. 

As one make the system more anisotropic, and moves towards the XYl phase the gap 
in the symmetric sector decreases. The critical values of J_l,ci become larger than J]_ for 
K > 0.76465. For such values of K ,the 2kF component of the random z-field is irrelevant. 
Therefore, for < J± < J_l,ci a gap is generated by interchain coupling in the antisym- 
metric mode but not in the symmetric mode due to the g ~ component of the random 
z field. In other words, an intermediate gapless XYl phase exists between the decoupled 
chain phase and the stable Haldane phase. Since the VBS order parameter is the correlation 
function (cos \/2(f)s{x) cos \/2(f)s{x')) in the limit | x — | (see appendix it is strictly zero 
in the phase that results from the breaking of the Haldane gap, due to the presence of the 
forward scattering disorder. This situation is shown in Figure ^b. If i^' > 1, is irrelevant 
in the symmetric sector and one recovers a direct transition between the decoupled chains 
and the stable XYl phase, with algebraic correlations, as shown on Figure |^c. 

For random exchange, it has been shown in Sec. |I V A 2| that the g ~ component of 
disorder is irrelevant in the RG sense. For values of K in the interval [0.76465, 1] the q ~ 2kf 
part of the disorder is irrelevant, so random exchange which does not give rise to a relevant 
g ~ term will continue to give the transitions of Figure ^a. For K > 1, random exchange 
disorder gives the transitions of Figure ^c. 

It is noteworthy that the resulting phase diagrams for disorder larger than the gaps in 
the ladder system is quite different from the one where no coupling of the antisymmetric 
mode to the disorder could occur such as for a true spin one chain. In that case increasing 
the disorder would lead to a disordered Haldane phase that could retain some degree of VBS 
order for the random exchange caseSil. 



V. RANDOM PERTURBATION BREAKING THE U(l) SYMMETRY 

From a theoretical point of view it is also interesting to consider randomness that break 
the U{1) rotational symmetry of the XXZ model. Indeed for the single chain casei, such 
disorder was proved to be very efficient in destroying the quasi long range order in the 
system. Given the remarkable stability of the gapless phase of the ladder system to the 
type of disorder examined in section it is interesting to check whether the same property 
still occurs. In this section, the two most common types of disorder breaking the U{1) 
symmetry are considered, namely a random field confined to the XY plane and random 
planar anisotropy. 



A. Coupling to disorder 
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1. Random field in the XY plane 



The two chain system with a random magnetic field in the XY plane is considered. The 
Hamiltonian is H = impure + Hxyf, Hxyf being given by 



(37) 



With hlph'^j^g = DSij5p^g5a,b ia,b = x,y p, q = 1,2). 

To bosonize this expression, it is convenient to rewrite (|37D in terms of Sf^^, (a = 1, 2) 
and introduce hf^ = h^^ ± th^a- This gives 



(38) 



For chain 1, upon bosonization the coupling with the random XY field can be rewritten sM: 

dx 



XYF,1 



Tia 



/il'^~°)+e-^^^(^)cos20i(a;) + /^f ^)+e-^i(-)l + H. c 



(39) 



A similar expression holds for chain 2. From power counting, the most relevant terms in 
the bosonized Hamiltonian are the 2kp ones. One could be tempted to keep only the 2kF 
terms and simply drop the g ~ ones. However, the g ~ part of the coupling to disorder 
has a bosonized form e^*^ cos(20) as can be seen from (^) and generates in second order in 
perturbation theory a term / dxhlf^^{x) cos{2(j)i){x) i.e an effective random magnetic field 

parallel to the z axis (see (§)) with hl^ -^{x) oc h^i^^^'^hf'''^ + h. c. That term is always 
more relevant than the g ~ term, and can be relevant even when the 2kp term is irrelevant. 
It is easily seen that no other relevant terms are generated so that one has to keep just the 
generated random z field and drop the g ~ term from the Hamiltonian. This situation 
is typical of random perturbations that break the rotation symmetry around the z axi^. 
The generated random 2;-field transforms the random field restricted to the XY plane into 
a more general anisotropically distributed random field. As a byproduct, the effect of an 
isotropic random field can also be easily obtained. It is obvious that the phase boundaries 
induced by the random XY field and the isotropic random field will be identical. However, 
the correlation length will not show the same dependence on disorder strength in the phases 
the properties of which are dominated by the z component of the random field. 

From the preceding discussion of the generated terms, the bosonized Hamiltonian con- 



taining only the most relevant terms is H 
and 



bosonized 

with ifpure given by (|TUp 



XYF, bosonized 



dx 



^ 1 ' e ^ + ^ 1 ' e V2 + H. c 

v27ra v27ra 

+ / dx [/l^ff. i(x) COs(V2(0a + 0s)) + /ieff. 2(2^) COs(v^(0a - 0s)) 



(40) 



where /iefr.,p(a;)/iefr.,p'l 



X' 



CD'^5{x — x')6p^pf The second line of the above equation is the 



random z-field already studied in Section the first line is the random planar field itself. 
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One can expect the phases showing strong XY fluctuations to be more strongly affected 
by the random XY fleld than the phases having dominant antiferromagnetic fluctuations. 
Indeed the latter ones only couple to the generated random 2;-fleld, or in other words their 
coupling to disorder involves intermediate states with energies above the gaps. Technically, 
this means that it is necessary to separate the case in which (pa is massive (i. e. dominant 
antiferromagnetic fluctuations) and the case in which 9a is massive (i. e. dominant XY 
fluctuations). Obviously, the singlet phase and the Haldane gap phase should be insensitive 
to small random flelds since these phases are completely gapped and have vanishing spin-spin 
correlations. 



2. Random planar anisotropy 

Another form of perturbation breaking the rotational symmetry around the z axis is 
when 6J = Jx — Jy fluctuates randomly from site to site with a mean value zero. The 
relevant Hamiltonian is impure + -f^RA in which 

hka = y. - + ^'^(5^:2^^1,2 - si,si,^,)\ (41) 



With 6Jf6Jj = D6ij6p^q {p,q = 1,2). As in the case of the random XY flelds, contractions 
of the g ~ and q 2kp terms generate relevant termsi. The bosonized form of the g ~ 
part is : 

H'ra = J dx^ii.^) cos{V2{ea + es)){x) + J dx^2{x) cos(v^(0, - 0J)(x) (42) 

with 7p(x = na) oc 6J^, and the bosonized form of the 2kp part is : 

^RA = / dx^l'^-ix) cos(y2(^a + es))ix) cos(V2(0, + 0,))(x) 

+ / rfa;72'^(a;) cos(V2(0, - ea))ix) cos(v^(0, - 0,))(x) (43) 

An effective random z fleld is generated by contraction of the g ~ with the 2kF part as in 
the case of the random fleld in the XY plane. A Its expression is 



-^Etenerated ' dx 



hlix) cos{V2{(Pa + (ps)) + hl{x) co^{V2{<Pa ~ 4>s))\ (44) 



where hp{x)hp,{x') = CD'^5{x — x')5p^pi This random z fleld ( PD is always more rele- 
vant than the 2k p part. Therefore only (^) and (^) have to be kept in the bosonized 
Hamiltonian. 



B. Physical Properties 

Having worked out the bosonized representations of the perturbations that break the 
f/(l) symmetry, it is possible to analyze their effects on the phases of the system. 
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1. XYl phase 



This phase is expected to be the most unstable in the presence of a random XY field or a 
random anisotropy. Since the XYl phase is not affected by a random z-field (see Section |rv|), 
the corresponding term can be safely dropped. 6a is replaced by its average value {6a) in 
(HO) giving the following simplified expression of the coupling to the random XY field: 



-f^CYF, eff. ^ I dx 



eeff.(a;)e +H. c. 



(45) 



Where ic?i.[x) = (^i(x) + G(x)) for J_l < and ^es.{^) = "^(^li^) ~ ^2{x)) for J± > (see 
tables 10). 

In the case of the random anisotropy, the effective coupling in the XYl phase has the 
form: 

^RA = / 7es.ix) smiV26six))dx (46) 

for J± < and: 

Hra = J lesXx) cos{V26s{x))dx (47) 

for J± > 0. 0s being massless, the RG equation for the random XY field is derived from 
(^) in the form 

(48) 



dl V 4:Ks 

Equation ( P5| ) implies that the random XY field is relevant for Kg > 1/12. Since the XYl 
phase only exists at i^'s > 1 (see tables |,|D this phase is completely suppressed by an 
arbitrarily weak random magnetic field in the XY plane. At small disorder the correlation 
length of the disordered XYl phase is 

k. field ~ (1/D)^5^ (49) 

Similarly for the case of planar anisotropy, the RG equation for the disorder is 

dD f I \ 

The disorder is thus relevant for Ks > 1/3 and the XYl phase is also unstable to weak 
random anisotropy. The correlation length is in that case 

k anisotropy ~ (51) 

For the isotropic point l^, anisotropy ~ 1/-D. 



20 



2. Singlet and Haldane gap phases 



In these two phases, the couphng with disorder has the same simphfied expression (H3-H^ 
and (1^) than in the XYl phase due to the identical structure of the gaps in the antisym- 
metric modes. Moreover, it is known from section |^ that a random z-field has no effect. 

However, the presence of a mass term in the symmetric part of the Hamiltonian leads to 
the suppression of all disorder correlation functions. The gapped phase is thus stable as for 
the case of perturbations preserving the rotational symmetry. 



3. The XY2 phase 

This phase has weaker XY fluctuations than the XYl phase and should therefore be less 
affected by the random XY field or the random anisotropy . Nevertheless, it has subdominant 
antiferromagnetic fluctuations and can be disordered by coupling to the generated random 
z field. 

Indeed it is known that random z fields suppress the XY2 phase if i^'s < 3 (see Section 
fVp. For Ks> ?) on the other hand, only the random XY field or the random anisotropy can 
possibly disorder the XY2 phase. In the XY2 phase, is massive, and this gives apparently 
zero coupling to disorder when simplifying (^Dp. In fact, an effective coupling of the random 
XY field to (f)s can be derived through second order perturbation theory along the lines of 
Ref. P3|J55|. The calculations are straightforward and lead to the following effective coupling 



to the random XY field 



HxYF, eff. = I dx 



U{x)e'^'' c.l (52) 



where C,ea.{^)^ea.{^') ^ D'^6{x — x'), and prefactors coming from mean values of the (pa 
fields have been omitted. A similar calculation for the random anisotropy case leads to the 
following effective coupling: 

^^RA, eff. = J dx'jes.ix) cos( V86',(x)) (53) 



With 7eff.(a;)7eff.(a;') oc D^6{x — x'). 

From power counting, the RG equation for the random field is 

dDxY /3 1 



dl \2 2K, 



DxY (54) 



implying that the 2kF part of the random XY field is relevant as soon as Kg > 1/3. The 
correlation length in the disordered phase is 

^random field ~ (V Dxy) ^^^'^ (55) 

For random anisotropy, the RG equation is: 
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'l-l)D.. (56) 



dl \2 K, 

Since the random planar anisotropy is relevant for Ks> the XY2 phase does not survive 
the presence of a small random planar anisotropy also for Ks> ?>. the correlation length of 

2 if 3 

the resulting disordered phase is then given by /random anisotropy ~ (l/-Di?A)^^=-*- 



4- The Ising antiferromagnet 
Since the antiferromagnetic phase is unstable with respect to a weak random z-field (see 



section ^Vj) , it will be completely suppressed here by the generated random z field both in the 
case of random anisotropy and random field in the XY plane. The corresponding disordered 
phase is a random antiferromagnet of correlation length / ~ l/D^ . This correlation length 
is just the one of the disordered XY2 phase "frozen" at K = 1. 
The full phase diagram is given on figure M. 



C. Physical discussion 

1. comparison with the one chain system 

In order to compare the two XXZ spin chain system with a single XXZ chain, the 
discussion has to be restricted to the antiferromagnetic, singlet and XYl phases since the 
XY2 phase does not obtain in a two XXZ chain system. The single chain system has a XY 
phase for < Jxy and an antiferromagnetic phase for > Jxy- The XY phase of the single 
chain is suppressed by a random field in the XY plane and the correlation length i^ lict. ~ 

2K 

[1/ Dxy)^'^~^ ■ The XYl phase in the two chain system is also suppressed, but : InR^ch. ^ 
/RF.ich. (see (^9|)). Similarly, it is also unstable in the presence of random anisotropy^!, with 

a correlation length /ra,i ch. ~ i^/ Dra)'''^"'^ in the disordered phase. Again this length 
is much larger than its ladder counterpart /RA,2ch. ^ /RA,ich. (see (|5T1) ). Thus, interchain 
coupling makes an XXZ ladder systems much more sensitive to random perturbations that 
break rotational symmetry around the z axis than a single chain. In both case the effect of 
disorder is considerably enhanced since the exponents of the correlation length are changed. 
By contrast, for perturbations that preserve the rotation symmetry around the z axis the 
one chain system was more disordered than its two chain counterpart. This seem to indicate 
that the XY phase in the ladder system is in effect much more anisotropic than its single 
chain counterpart. 

As for perturbations respecting the rotation symmetry, effects of disorder in the anti- 
ferromagnetic phase are quite similar for a single chain and for the ladder. Indeed for a 
single chain, contrarily to the results of Ref. ^ , the generated random z field suppresses the 
antiferromagnet phase, both in the case of a random field in the XY plane and of random 
anisotropy. The same effect also holds in the two chain system. In that case, the correlation 
lengths show the same dependence with disorder strength up to prefactors. Contrarily to 
the one chain case, there is no phase boundary between the disordered XY phase and the 
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random antiferromagnet phase due to the existence of the Haldane gap (or the singlet) phase 
which is not suppressed by a weak disorder. 

The XY2 phase is also suppressed. For a random field, however, the correlation length 
in the disordered phase is much longer than the one of the disordered XY phase of a single 
chain, /corr. ~ (l/D)<3^-i. In the presence of random anisotropy, the XY2 phase is also 
suppressed. However, the XY2 phase gives rise to two different disordered phases. The 
first one corresponds to suppression of the quasi long range order by the random anisotropy 
per se, the other one to suppression of QLRO by the effective random z field. A crossover 
occurs between these two phases at the point where they have identical correlation lengths. 
The correlation length of the phase pinned on the random z field being given by /corr. ~ 

2 

{1/ Dz)^-^^ where measures the strength of the effective random magnetic field. Since 
Dz ~ the crossover should occur for Kg = 2. for Kg < 2, the properties are dominated 
by the generated random z field and for Kg > 2 by the random anisotropy. Two different 
disordered phases already existed in the XXZ chain with random planar anisotropyl^. Such 
an effect in the two chain system is specific of the XY2 phase and does not exists in the 
XYl phase. The correlation length of the disordered phase induced by a random anisotropy 
in a XY2 phase is larger than the correlation length of its one chain counterpart /corr.,i ch. ~ 
3ifi^. Paradoxically, although the XY2 phase does not exist in an XXZ two chain 
system with weak antiferromagnetic coupling, its properties regarding disorder are much 
more similar to the one of the XY phase of the single XXZ chain than the ones of the XYl 
phase. Such difference is explained by the fact that in the XYl phase the 2fci? fluctuations 
in are suppressed by gap formation leading to an increased robustness with respect to 
perturbations that respect rotational symmetry around the z axis and by the locking of 
the spins that sit on the same rung of the ladder leading to an increased sensitivity to 
perturbations that break rotation symmetry around the z axis. Also, it can be shown that 
in the XY2 phase, the spins on each rung of the ladder can only take identical values (for 
J± < 0) or opposite values (for J± > 0) leading to an effective spin 1/2 chain. 



2. Effect of a uniform magnetic field 
As seen in section [IV C 2|, a uniform magnetic field can inhibit Haldane gap formation. 



In the presence of a random XY field or random anisotropy, this lead to the formation 
of a disordered XY phase, with a correlation length given by Eq. ( |^ ) for a random XY 
magnetic field and by Eq. (|5TD for a random planar anisotropy. It can be verified that the 
correlation length induced by the effective random z field is much longer than the correlation 
length induced by the random XY field for Kg > ~ 0.106 and much longer than the 

correlation length induced by the random XY anisotropy for Kg > ~ 0.3802. In other 

words, the induced random z field is a very weak perturbation that takes over only when the 
random field in the XY plane or the random anisotropy is nearly irrelevant. Therefore, in 
the presence of an applied uniform magnetic field, a disordered phase replaces the Haldane 
gap phase. The resulting phase diagram can be found on figure |^ for the random magnetic 
field in the XY plane and on figure |^ for the random planar anisotropy. 
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VI. CONCLUSION 



In this paper we have investigated the effects of disorder on an anisotropic two leg spin 
ladder. 

For the pure system, we have computed the phase diagram in details. Both for ferro- 
magnetic and antiferromagnetic interchain coupling the system can exhibit four different 
phases: an antiferromagnctically ordered state, a gapped singlet phase (Haldane phase) and 
two massless XY phases. One of the XY phases (XYl) is a close analogous to the XY phase 
of a single chain. The other one (XY2), that was not discussed in previous work on the 2 
coupled chain system, is the analogous to the one occurring for a spin one chain with on site 
anisotropy. That phase has a tendency to have the spins parallel to the z axis. For weak 
interchain coupling only the AF, Haldane and XYl phases can be realized. An interesting 
question is whether for intermediate coupling one can stabilize or not the XY2 phase. Of 
course such a phase can always be realized by using more complicated interchain couphngs. 

The disordered ladder shows remarkable features compared to a single spin chain. For 
perturbations respecting the XY symmetry such as a random 2;-ficld and random exchange, 
both the singlet phase and the massless XYl phases revealed to be totally insensitive to 
weak disorder. On the other hand the XY2 phase is extremely sensitive to disorder due to 
the presence of strong antiferromagnetic fluctuations in the z direction. Similarly to the 
one chain case, the effect of randomness on the Ising antiferromagnet depends on whether 
the perturbation is invariant under ^ —S^ or not. In the former case it does not affect 
the Ising AF in the latter it suppresses it through an Imry Ma mechanism. While such a 
stability was to be expected from the gapped singlet phase, it is much more surprising in the 
massless XY phase and is reminiscent of the delocalization occurring in fermionic ladders. 
Note that here the effect is even stronger, since the delocalization transition only occurred 
for weakly attractive interactions close to the non-interacting point. This would correspond 
to a non-disordered phase for Jz < and a disordered phase for > in the spin language, 
whereas the spin system is stable with respect to disorder for Jz < J. In the presence of a 
finite magnetic field this stability to disorder subsides, even if the singlet phase is destroyed. 
The isotropic point K — 1/2 becomes now sensitive to disorder, but for moderate XY 
anisotropy K > 3/4 disorder has again no effect. This problem offers interesting connections 
to disordered bosonic ladders or coupled vortex planes. This extreme stability to disorder 
for the spin ladder is to be strongly contrasted with the one for a single chain where most 
of the phase diagram is destroyed by infinitesimal disorder and only extremely anisotropic 
XY chains K > 3/2 can resist. 

Knowledge of the phase diagram for weak disorder allows to make reasonable guess 
on the behavior of the system for stronger disorder or when the gaps are reduced (e.g. 
by diminishing the interchain couplings). Close to the isotropic point K ~ 1/2, upon 
increasing the interchain coupling, there should be a single transition between a regime of 
decoupled (and thus disordered) spin 1/2 chains towards a stable singlet phase ladder for 

~ £)(2-i/2/sr)/(3-2K) However by making the system more XY anisotropic the transition 
should occur in two step. Decoupled spin chains would have a transition towards coupled 
chains, where the antisymmetric mode is gapped, but the symmetric mode is not gapped, 
leading to an XYl phase. Upon increasing the interchain coupling further the symmetric 
mode gaps giving back the singlet phase. Thus for anisotropic ladders with 0.76465 < K < 
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1 an intermediate non-disordered XYl phase should appear for intermediate interchain 
couphngs. For large anisotropics the singlet phase does not exist any more and one recovers 
a direct transition between decoupled chains and the stable XYl phase of the ladder. Since 
this transition are obtained by the crude comparison of the various correlation length in the 
system it would be interesting to confirm such a phase diagram by numerical simulations. 
Subtle effects might indeed occur for disorder such as the random exchange where the true 
correlation length of the systemf compared to the one given by the RG around the Gaussian 
fixed point) is found to diverges. 

The extreme stability of the ladder to XY symmetric randomness made it worth to 
investigate perturbations breaking this symmetry as well. Remarkably the behavior is here 
inverted. Randomness breaking the rotational symmetry around the z axis, suppresses the 
two XY phases and the Ising AF phase as was the case for one chain. However the disordered 
XYl phase has a much shorter correlation length than the disordered XY phase of the one 
chain system in the presence of the same perturbations whereas the XY2 phase has a much 
longer correlation length than the XY phase of the one chain system. Such an effect is 
again due to the fact that the ladder system with a given Jz/J 7^ 1 is in effect much more 
anisotropic than its one chain counterpart with the same Jz/ J parameter. The XY2 phase 
in the presence of planar anisotropy gives rise to two different disordered phases similarly 
to the XY phase of the one chain system but at odds with the XYl phase that gives rise 
to only one phase. Such result is at first sight paradoxical since the XYl appears as the 
natural continuation of the XY phase of the one chain XXZ system, whereas the XY2 phase 
is likely to exist only in more complicated models or at larger coupling. 

In the presence of a strong enough uniform magnetic field, the Haldane gap phase is 
suppressed, allowing for the observation of a crossover from a phase dominated by the 
random field in the XY plane or the random planar anisotropy to a phase dominated by a 
generated random field parallel to the z axis. This crossover occurs in the vicinity of the 
transition to the disordered classical antiferromagnet and may therefore be difficult to probe 
in experiments or numerical simulations. 

Clearly the disordered ladder system presents an extremely rich behavior and the renor- 
malization study presented here can only be a first step towards its understanding. This is 
specially true for the crossover to strong disorder, which would be interesting to investigate 
through numerical simulations and non perturbative methods. 
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APPENDIX A: SINE-GORDON HAMILTONIANS 

The general sine-Gordon Hamiltonian has the form: 
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27r 



u 
K 



A dx cos(/50) 



(Al) 



To get some understanding of let us first put A = 0. The simphfied Hamiltonian is 



H= I — 

2-K 



(A2) 



Computing the correlation functions at OK for (|A2|) in Matsubara time gives 

^rp ^ine{x,T)^-m8{0,0)^ ^ j + {utY 



\2\ — 



\2\ -2K 



(A3) 
(A4) 



Thus K controls the power law decay of the correlation functions i. e. the scaling dimensions 



of the operators. More precisely 



imply that e*"'^ has dimension and that e 



has dimension u can be interpreted as the velocity of the excitations. To compute the 
spin-spin correlation functions, one uses the expressions (^. 

For A 7^ 0, correlation functions cannot be obtained exactly anymore. However, the sine- 
Gordon Hamiltonian can still be studied using renormalization groupfRG) technique^lS. 
The flow equations for K and A are of the Kosterlitz-Thouless forniEal^. From ([A^), A has 
scaling dimension 2 — j3'^K/A. A small A is thus irrelevant for K > Kc = 8//3^. 

When A is irrelevant, K flows to a fixed point value K* and correlation functions keep 
their power law character up to logarithms!^ with the bare value of K replaced by K*. On 
the other hand if A is relevant, acquires an expectation value that minimizes the ground 
state energy and a gap is formed. It can then be shownS that there {f{4>)) ~ /((0)) and 
,zae{x,T)^-taem^ _ cxp( 



that {T^e'' 



where ^ is a correlation length. These results 



are used extensively in the paper. 



APPENDIX B: CALCULATION OF THE VBS ORDER PARAMETER IN THE 
PHASE IDENTIFIED AS A HALDANE GAP PHASE 

In this appendix the identification of sector II of table | with a Haldane gap phase is 
made more precise. It is a well known fact that a Haldane gap phase has a hidden topological 
long range orderil'il . The order parameter that measures the hidden topological orderel is 
known as the Valence Bond Solid (VBS) order parameter. 

The VBS order parameter C is a nonlocal order parameter defined as 




In the Haldane gap phase, all the spin-spin correlation functions decay exponentially but 
C 7^ 0. Physically, the fact that the VBS order parameter is non zero indicate that if all the 
sites with 5*^ = are removed from a spin-1 antiferromagnetic chain 
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the remaining ("squeezed") chain has antiferromagnetic order. 

In the following, it is shown using bosonization that the mean value of the VBS order 
parameter is non zero in the phase we identified as a Haldane gap phase. A related problem 
is the dimerized spin chainll, where a bosonization transformation permits to show explicitly 
the existence of hidden long range order. 

In the first place, a bosonized expression of the string operator exp{ni J2i<n<jiSi + S2)) 
has to be obtained. A naive bosonization would lead to a string operator of the form: 

{ f^iA.c\ P'^ A cos(y20,) cos(y20a) 
exp(2V2(0s(ja) - 04«a) + / dxe )) (B2) 

Jia Ha 

That expression has very bad features. First, in the antiferromagnetic phase, the oscillating 
term is non zero and creates a complicated variable phase, that is extremely tricky to handle 
because of the unknown lattice renormalizations . In particular it is impossible with that 
expression to get the mean value of the VBS order parameter. Moreover, in the Haldane 
phase the mean value of the oscillating term is zero but if it is dropped of the calculation of 
the VBS order parameter , an unphysical zero VBS order parameter results. However, these 
incorrect results are only due to an inappropriate boson representation of the string order 
parameter. The good representation is obtained using the identity: exp(«7r(S'^ + 5*2)) = 
— exp(z7r(5'f — 5*1)) .This operator identity comes from the fact that = ±1/2 so that the 
string order parameter can now be rewritten 

n exp(m(5i^ + S^)) = i-y—' exp( ^ zniS^ - S^)) (B3) 

i<n<j i<n<j 

The new form of the string order parameter can be straightforwardly bosonized in the form 
(—)*"■' e*^'^'^"*^*"^"''^''*^-'"^^ This one has the correct sign alternation in the antiferromagnetic 
phase. Using the bosonized expression of 81+82 = 8^ 

5^(3,) _ -V^^x0. ^ ^^^x/a 2 C0S( V^0,) COs( ^20 J ^^^^ 

The VBS order parameter is obtained in the form 

C= lim (cos(y20,(x))cos(y20,(y))) (B5) 

Since (ps is long range ordered, a non zero VBS order parameter is obtained in the phase 
identified as a Haldane gap phase in section [HI A| . More precisely, C ~ ((cos(v^0s)))^. An 



alternative derivation can be found in Ref. ITS 



APPENDIX C: DERIVATION OF THE BOSONIZED COUPLING FOR A 

RANDOM Z EXCHANGE 

For the sake of simplicity, the derivation is made for one chain. 
It is done in two steps. First, one goes to the continuum limit 

-f^randomJz= / dxJ^{x)8''{x)8''{x + a) (CI) 
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with Jz{x) = aJn,x = na,S^{x) = Sn/a. 

Then, one fermionizes using the Jordan Wigner transformation ^ 

After straightforward calculations the 2k p part of S^{x)S^{x + a) is obtained in the form 

- {^PlijR + + ^l^b){x + a) + {^Ul + ^l^B)ix)iijUR + ^\i\)l){x + a) 

(C2) 

Written as above, this is a purely formal expression that contains hidden infinities. For it to 
make sense, it must be normal ordered . note that in the preceding cases, only two fermion 
operators had to be normal ordered . Thus, normal ordering gave only some constants that 
could be safely dropped from the Hamiltonian. Here, a product of 4 fermion operators has 
to be normal ordered. Thus, contractions can leave us with generated 2 particle operators 
that cannot be obtained in a naive bosonization procedure. 

Technically, Wick theorem says that any product of operators can be written as a sum of 
normal ordered products with pairings^. The two non zero pairings are the following 



2t\{x — x' 
1tx{x — X 



(C3) 



Let us pick up the first term with a non-zero contraction 

+ a)i)L{x + a) = : V^|j(x)V'/j(a;)V'|j(x + a)i)L{x + a) : 

+ : i)^j^{x) i)f.{x)i)'^^{x + a) i)L{x + a) : 

= : iIj'^j^{x)iIjr{x)iIj'^r{x + a)?/'L(a; + a) : +— ^ : + a) 

i. e. upon normal ordering, the short distance expansions of operator products generates a 
2 fermion term that cannot be obtained doing a naive bosonization. 

Looking carefully at expression one can see that it contains 4 generated 2 fermion 

terms plus the fully normal ordered 4 fermion terms. 

Collecting together the 4 generated 2 fermion terms, one gets 

^: + a) : - : iJL{x)ipR{x + a) : + : ip\i{x)tpL{x + a) : - : ^/'l(a;)^R(a; + a) 



2«7ra 



(C5) 

Upon bosonization, this gives 

J_ X -L (e-* - e-'A = (C6) 

Thus, one obtains 

(ia;— 7 — TT— sin 2(\){x) + less relevant terms . . . (C7) 
vra r 
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The normal ordered 4 fermion terms reduces to a term {dx(f))'^ sin 2(f). Power counting 
implies that that term has dimension 2 + K and is thus irrelevant irrespective of the value 
of K. 

It is possible to check that (|C^) is the correct expression by the following argument: 
S^S^j^i is invariant under a rotation around the z axis, so that it cannot depend on 9 and since 
it must have the same scaling dimension at the isotropic point as S^Sfj^^, its bosonized form 
must be either cos 20 or sin 20. Since —>■ —S^ corresponds in the bosonization language 
to —> I — the bosonized form of SfSfj^i must be invariant under such transformation. 
This rules out the cos 20 and leads to ( P7| ) as the only possible expression for the bosonized 
form of S■S^_^_l. However, it is comforting to be able to derive directly the expression (|C^), 
since it proves that (|C7|) is not an ad hoc expression that is chosen arbitrarily in order to 
obtain physically sound results from an inappropriate technique. 
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TABLES 



TABLE L the 4 sectors of the pure 2 spin chains model with J± < 





I 


II 


III 


IV 


Ks 

Ka 


< 1 
< 1/2 


< 1 
> 1/2 


> 1 
< 1/2 


> 1 
> 1/2 


4>s 

Oa,<t>a 

phase 
Order Parameter 


(0s) = 
(0a) = 

Ising AF 
cos(a/2^s) cos(\/2^a) 


(0s) = 

{ea) = 

Haldane gap 
? 


massless 

(0a) = 

XY2 


massless 

{ea) = 

XYl 
cos(^) 


TABLE IL the 4 sectors of the pure 2 spin chains model with J± > 




I 


II 


III 


IV 


Ks 

Ka 


< 1 

< 1/2 


< 1 

> 1/2 


> 1 

< 1/2 


> 1 
> 1/2 


<t>s 

phase 
Order Parameter 


= ^ 

(0a) = ^ 

Ising AF 
sin(A/2^s) sin(v/20a) 


^<t>s) = ^ 
{0a) = ^ 

singlet 
VBS 


massless 

(0a) = ^ 

XY2 
QiV2es 


massless 

{ea) = ^ 

XYl 
sin(-^) 
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FIGURES 
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FIG. 1. The phase diagram of the pure 2 chain XXZ model in terms of and Kg XYl and 
XY2 are gapless phases (see text). AF contains antiferromagnetic quasi- long range order. The 
singlet (antiferromagnetic interchain coupling) or Haldane (ferro. interchain coupling) have a gap 
to all excitations. The dotted line represent the weak interchain coupling case, when the intrachain 
anisotropy is varied. The isotropic point is Kg = Ka = 1/2. The dashed line corresponds to two 
coupled XXZ chains. 
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FIG. 2. The phase diagram of the 2 chain XXZ model with a random z field in terms of and 
Kg. Lines indicate the parts of the phase diagram where disorder is relevant. The singlet phase 
is stable due to the presence of the gap. Quite surprisingly the massless XYl phase is now also 
totally insensitive to weak disorder. The dashed line is for two coupled XXZ chains. 
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FIG. 3. Phase diagram for fixed disorder as a function of the strength of the interchain coupHng 
Jl- (a) For K = 1/2 the isotropic point. There is a direct transition between decoupled disordered 

spin 1/2 chains and a stable Haldane or singlet phase, (b) For K > 0.76465 a disordered XYl 
phase exists between the decoupled chains phase and the stable Haldane phase. This phase shows 
no topological order, (c) For K > 1 A direct transition between decoupled chains and the stable 
XYl phase takes place. 
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FIG. 4. The phase diagram of the 2 chain XXZ model with a random planar exchange or a 
random z exchange in terms of Kq, and Kg. Lines indicated the parts of the phase diagram where 
disorder is relevant. Here again the XYl phase remains unaffected by weak disorder. The dashed 
line represents coupled XXZ chains. 
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FIG. 5. The phase diagram of the two chain XXZ model with a random z magnetic field under 
a uniform magnetic field parallel to the z axis. The uniform magnetic field is strong enough to 
inhibit Haldane gap or singlet gap formation. This results in the apparition of a disordered phase. 
The dashed line corresponds to coupled XXZ chains. 
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FIG. 7. The phase diagram of the two coupled spin chains in the presence of a random field 
in the XY plane and a uniform field parallel to the z axis. The dashed line corresponds to two 
coupled XXZ chains. The uniform field inhibits Haldane gap formation. 
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FIG. 8. The phase diagram of the two coupled spin chains in the presence of a random pla- 
nar anisotropy and a uniform magnetic field parallel to the z axis preventing Haldane gap phase 
formation. The dashed line corresponds to two XXZ coupled chains. 
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